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Deformation quantization produces families of mathematically equivalent 
quantization procedures from which one must select the physically mean¬ 
ingful ones. As a selection principle we propose that the procedure must 
allow enough ‘observable’ energy distributions, i.e., ones for which no 
pure quantum state will appear with negative probability and must fur¬ 
ther have the property that for these the uncertainty in the probability 
distribution of the quantum states must not exceed that of the original 
distribution. For the simple harmonic oscillator we show that this allows 
only the classic Groenewold-Moyal (skew-symmetric) form. 

The idea of negative probabilities is not new, going at least as far back as 
Wigner and perhaps even into the 19th century. It has been likened popularly 
to observers on the sidelines of a soccer game seeing (non-negative) probability 
distributions of the coordinates of the ball, none of which may seem unusual, 
but from which they deduce that the probability distribution of the ball over 
the interior of the entire playing field has points where it is negative. Suppose, 
however, that we have a system for which we know some distribution of ener¬ 
gies. Deformation quantization generally produces a family of ‘cohomologically 
equivalent’ quantizations of the system, each of which together with the en¬ 
ergy distribution assigns a probability, possibly negative, to each pure quantum 
state. An energy distribution will be called observable with respect to a given 
quantization procedure if each of these probabilities is in fact non-negative. We 
will say that there are enough observable distributions if every distribution is 
in the closure of the linear space spanned by the observable ones. 

With any quantization procedure one can associate to the original distribu¬ 
tion of energies two measures of uncertainty (=standard deviation), that of the 
original distribution and that of the distribution of energies it produces in the 
various pure quantum states. The latter will be called the quantum uncertainty ; 
The selection principle proposed here is that amongst cohomologically equiv¬ 
alent quantizations only those are physically meaningful for which there exist 
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enough observable distributions, and where for each observable distribution the 
quantum uncertainty does not exceed that of the original distribution; briefly 
quantization must decrease uncertainty. 

In the case of the simple harmonic oscillator, examining the entire family of 
quantizations possible through deformation yields identities involving Laguerre 
polynomials which, aside from the present method of derivation, are generally 
not new. They show, however, that for each of these quantizations there is 
a natural infinite sequence of ‘basic’ observable distributions from which any 
pure state can be recovered as a linear combination. In particular there are 
always enough observable distributions, so for the simple harmonic oscillator 
this by itself is no restriction on deformation quantization. However, the only 
quantization that the inequality on uncertainties allows is the Groenewold- 
Moyal form (although the normal or anti-normal form may be meaningful when 
quantizing fields). To single out that form in the case of the simple harmonic 
oscillator it would be sufficient, as will be shown, to require that the quantum 
uncertainty of an observable energy distribution approach that of the original 
distribution as the energy tends to infinity; perhaps this alone would be sufficient 
in general. 

1 Some basic algebraic deformation theory 

The seminal paper in deformation quantization is that of Bayen, Flato, Frpnsdal, 
Lichnerowcz and Sternheimer [I] 1978, some essential ideas and results of which 
are simply reproduced here without further attribution. This approach to quan¬ 
tization has been exceedingly fruitful. Some of the subsequent developments are 
summarized in [2j , which is complemented by an extensive bibliography. (A use¬ 
ful recent introduction to the theory is the note of Hirshfeld and Henselder {2j.) 
We begin with a very brief review of algebraic deformation theory, introduced 
by the author in j3j. 

Let A be an algebra which here (and generally in any physical theory) will be 
assumed to be over the real or complex numbers, e.g., the algebra of functions on 
a phase space, but in a more general context could be over an arbitrary commu¬ 
tative, unital ring. A deformation of A is a new associative ‘star’ multiplication 
expressible as a formal power series 

a * b = ab + HC\{a , b) + H 2 C 2 {a, b) + .... (1) 

Here h is for the moment just a formal parameter and the C) are bilinear maps 
from A x A to A. We should like the star product to be defined on the same 
underlying vector space as A but the introduction of the formal parameter h 
generally makes it necessary to extend the coefficients. Frequently the extension 
is to power series in H and one views the functions C\, C 2 ,... as having been 
extended to be bilinear not only over M or C but also over these power series. 
Often one tacitly assumes this done and in favorable cases the series actually 
converge for sufficiently small values of the variable H. In deformation quanti¬ 
zation, however, this extension of coefficients is not appropriate because we are 
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forced to consider power series in 1 /H. The correct extension of coefficients in 
this case is to the field of Laurent series in 1/H, i.e. C[[l/fi]][/i], but this puts 
additional restrictions on the Ci in order for the star multiplication to be mean¬ 
ingful. It is sufficient that the star product 0 be locally finite , i.e., that for 
any a,b £ A only a finite number of the Cfia , b) be non-zero. (This is likewise 
frequently unmentioned but it is generally automatically satisfied when the Ci 
are bidifferential operators). The use of Laurent series in 1 /h rather than power 
series in h has, as we will see, a profound effect on the structure of the resulting 
algebra. 

If we set Co (a, 6) = ab , the original associative multiplication, then associa¬ 
tivity of the star multiplication is equivalent to the condition that 

[Ci(Cj(a,b),c) — Ci(a,Cj(b,c))] = 0, all n > 1; *, j > 0 

i-\-j=n 

for all a, 6, c in the algebra A. Transposing to the right side all terms with i or 
j equal to 0, this becomes 

E [Ci(Cj(a, b), c) - Ci(a, Cj(b, c))] = 

i+j=n\ i,j>0 (2) 

aC n (b,c) — C n (ab,c) + C n (a,bc) — C n (a,b)c, all n > 1. 

These are generally difficult conditions to meet. In the Hochschild cohomology 
theory, each Ci is a 2-cocliain of A with coefficients in itself, and the right side 
in m is the Hochschild coboundary of C n . For n = 1 the left side is zero, so the 
coboundary of C\ is zero, that is, C\ is a 2-cocycle. This (or more properly, its 
cohomology class) is often called the infinitesimal of the deformation. A basic 
problem, given an infinitesimal deformation, is to construct a deformation which 
has it for infinitesimal. For n = 2 the left side of 0 is something constructed 
from Ci which in fact will always be a 3-cocycle in the Hochschild theory, and 
the first requirement is that it be a coboundary, namely the coboundary of Ci. 
Having Ci and C%, the left side of @ with n = 3 will again be a cocycle which 
is required to be a coboundary, and so forth for all n. Unless we have some 
control of the Hochschild cohomology of A it is clear that the construction of 
deformations in such a step-by-step manner will not be easy. 

Fortunately, there is one case in which all the conditions for associativity are 
automatically satisfied. A derivation I? of A is a linear mapping of A into itself 
such that D(ab) = ( Da)b + a(Db). We have been careful with the order of the 
variables a and b here because the multiplication in the algebra A need not have 
been commutative, although that is the case for any usual algebra of functions. 
Ordinary differentiation in the algebra of infinitely differentiable functions on R 
is a derivation. Suppose that D' and D" are commuting derivations of A. Then 
the star multiplication defined by 

a*b = ab + HD'a D"b + D' 2 a D" 2 b + -D ,3 a D" 3 b + ... 
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will be associative, cf. 2! > something easily verified by direct computation. This 
(despite the consternation of mathematicians) is frequently written as 

a* b = a e hD D b. (3) 

In the star product of m the C n are then given by C n (a, b) = (1 /n\)D' n aD" n b. 
More generally, if D\ ,..., D' r , D'{, ..., D" are all mutually commuting deriva¬ 
tions then 

is again an associative multiplication. 

A basic example illustrating 10 is that where A = C[q,p],D' = d q ,D" = id p . 
Then q*p = qp+ih , while p*q = pq , so [q,p\* = q*p—p*q = ih. The last equation 
determines the structure of the deformed algebra up to isomorphism but not the 
deformation (which contains more information), as different deformations can 
give isomorphic algebras. In fact, referring to the basic equation 0 , suppose 
that T : A —» A is a linear map of the underlying vector space onto itself of 
the form Ta = a + firia + H 2 T 2 a + .... Defining a *' b = T~ 1 (Ta * Tb) and 
denoting by A*, A*' the algebras with these two multiplications, the map T is 
an isomorphism A*> —> A*, so the new multiplication is also associative. We 
say that the deformations given by * and *' are cohomologically equivalent (c- 
equivalent), the adjective emphasizing that despite the isomorphism there may 
be some physical differences between the results. In particular, consider 0 
and take T = T\ = e xhD D where D'D" is just the composite of the two 
derivations and A is an arbitrary constant. Then it is easy to check that the 
resulting deformation is given by 


a*\b = a e t>D'D"-\h{D'D"-D"D') b ( 4 ) 

In particular, for A = 1/2 we have 

a Nqp b = ae%( D D ~ D D ^ b. (5) 

In the context of quantum theory, the deformation given by the commuting 
derivations D', D" in © is generally called the ‘normal’ form, that obtained by 
interchanging D' and D" or by setting A = 1 in 0 is the ‘anti-normal’ form, 
and the skew symmetric form of © is the ‘Groenewold-Moyal’ form (GM) 1 
m d Historically, the fact that cohomologically equivalent deformations may 
not be physically equivalent immediately raised the problem of selecting from a 
family of c-equivalent deformations those which are physically meaningful. This 
was already addressed in pQ where one important reason cited for preferring the 
GM form is its greater symmetries. Another is homological: Every 2-cocycle of 
C [q,p\ with coefficients in itself can be written uniquely as a sum of a symmetric 
part and a skew part. Both parts are again cocycles, but the symmetric part 

1 While often attributed solely to Moyal, the basic idea is present earlier in the work of 

Groenewold, and anticipated even earlier in works of Wigner and Weyl. 
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is always a coboundary; there is, up to constant multiples, a unique skew 2- 
cocycle, and that cocycle is a biderivation, i.e., a derivation as a function of 
each argument. As mentioned, in the case of the simple harmonic oscillator our 
selection principle allows only the GM form, but like the foregoing principles 
(symmetry, cohomological uniqueness) our selection principle should apply to 
many other cases. (The normal and anti-normal forms are excluded in the case 
of a single oscillator but may be essential when one has infinitely many, as when 
quantizing a field, and one must normalize the lowest energy level of each to 
zero to avoid having that of the whole be infinite.) 

Consider now the choice of coefficients of a deformed algebra. If we do 
not have some specific information about the 2-cochains C 1 ,; in 0 other than 
that they give a deformation (or about the derivations D' and D" in 0 other 
than that they commute), then for coefficients one must take the power series 
ring C[[ft]], else the formulas will not be meaningful. In this classic approach 
some basic algebraic properties of A are preserved. In particular, if a non-zero 
element a of A is not a zero divisor in the original multiplication (i.e., if there 
is no b ^ 0 such that either ab = 0 or ba = 0) then a will not be a zero divisor 
in the deformed algebra. Similarly, if a was invertible then it will continue to 
be so. It follows that a deformation of an integral domain will continue to be 
an integral domain, and a deformation of a division ring (= skew field) will 
again be a division ring. A deformation of a unital algebra remains unital and 
the deformation will be, in fact, c-equivalent to one in which the original unit 
remains the unit. However, a deformation of a commutative algebra like C[q,p] 
need not remain commutative; this is the basis of quantization. 

As mentioned, in a classical algebraic deformation we generally hope that 
the power series which are encountered actually converge for sufficiently small 
values of the deformation parameter H , but for purely algebraic purposes this 
may not be necessary. Suppose now, however, that the deformation has the 
local finiteness property that for every a and b in the original (undeformed) 
algebra there is an N such that Ci{a,b) = 0 for all i > N. The deformed 
algebra will then already be defined over the polynomial ring C[h), and we can 
extend coefficients, if we wish, to the field of Laurent series in 1/h. This is 
the case, for example, with A = C[g,p] and D' = d q = d/dq,D" = id p in 
0 . We must now also be careful in the definition of cohomological equivalence 
to require that T also be locally finite, i.e., that for all a there is an N such 
that Ti(a) = 0 for i > N. That is certainly the case for the T which gives the 
equivalence between the normal and Groenewold-Moyal deformations of C[q,p\ 
(with D' = d q , D" = d p ). 

When, in the locally finite case, we extend coefficients to Laurent series 
in 1/H the structure of the deformed algebra may be very different from that 
obtained with power series as coefficients. What was before deformation an 
integral domain may acquire infinitely many orthogonal idempotents; as a result 
there may be no natural way to apply a contraction in the sense of Inonii and 
Wigner [Bj to recover the original algebra. While this is inherent in deformation 
quantization, it raises difficult purely algebraic questions about the structure of 
the deformed algebra. The same is true, of course, of all algebras obtained by 
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c-equivalent deformations since they are all algebraically isomorphic. 


2 Deformation quantization of the simple har¬ 
monic oscillator 


The foundational paper pQ showed, in particular, that quantization of the simple 
harmonic oscillator could be viewed as an exercise in the deformation of the 
polynomial ring C[q,p], where now q and p are viewed as the position and 
momentum coordinates on the phase space M 2 . The Hamiltonian function for 
the simple harmonic oscillator is 


H(q,p) 


P 

2 m 


-Q ■ 


( 6 ) 


In the classic approach to quantization one substitutes for q and p operators Q 
and P which satisfy the fundamental commutation relation [Q, P ] = ih. Gener¬ 
ally this involves some ambiguity, for the Hamiltonian may contain monomials 
of positive degree simultaneously in p and q, but that is not a problem here. We 
can take, e.g., Q =multiplication by g, P = — ihd q and with this seek solutions 
to the Schrodinger equation 

ihil’ = H ip. 

This in effect chooses a specific representation of the Weyl algebra C{q,p}/(qp — 
pq — ih) and with this choice ip is viewed as a function of q and t. (The present H 
is time independent.) Mathematically this does not yet introduce any quantiza¬ 
tion; the latter is forced by the physical requirement that ip be square integrable 
with absolute value tend to zero at ±oo. By contrast, the deformation approach 
chooses a deformation which gives rise to the Weyl algebra, e.g. that in ® (i.e., 
such that the commutator of q and p is essentially their Poisson bracket) and 
rewrites the Schrodinger equation in the form 


ihip = H *ip. 


(7) 


This will be called the “deformation-Schrodinger” or d-Schrodinger equation. 

Were H a matrix operating on a vector ip, the solution would be e~ lHt ^ h ip{ 0). 
Here, bearing in mind that H is now an element of a non-commutative alge¬ 
bra with multiplication * one must still compute the exponential exp*(— iHt/h) 
where exp* indicates that the exponential must be computed using the de¬ 
formed multiplication. The problem is to express the result, which is an el¬ 
ement of the underlying vector space of the original (undeformed) algebra of 
functions on phase space, without reference to the deformed multiplication. 
There are now different possible choices for the deformed multiplication * but 
exp*(—iHt/H)ip(0) will always be a solution to the d-Schrddinger equation. 

In deformation quantization generally, one knows from pQ that 

exp *(-itH/h) = ^ 7 t e e ~ ltE l h (8) 

E 
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where the sum in the Fourier-Dirichlet series on the right is over the allowable 
energy levels E and the tte are functions on the phase space which are orthogonal 
idempotents in the * multiplication whose sum is 1. One has H * tte = Ette- 
(Note that we have tacitly extended coefficients to Laurent series in 1/h and this 
has introduced zero-divisors into the new algebra.) Further, these idempotent 
functions when integrated over all of phase space will yield a common constant 
which, in the case of simple harmonic motion whose phase space is the q,p 
plane, is 2ttH. We should like to interpret the tte, which are functions on phase 
space, as giving a probability distribution there but in general they may take 
on negative values (depending on the deformation chosen). Although negative 
probabilities may not be directly observable, we shall see that at least in the 
case of simple harmonic motion a reasonable interpretation as probabilities may 
still be possible. 

Following an idea often credited to Dirac, it is convenient to transform the 
Hamiltonian © into “holomorphic coordinates” by setting 


a = 


TiUJ p 

tt(9 + *-). 

2 mw 



■ p \ 

i - . 

mw 


With this one has 

H = uaa. 

The simplest deformation of C[a,a] one can now choose is the normal form 
defined by setting 

/ *n 9 = fe nSa °*g. 

One then has [ 0 , 0 ]*^ = h , which is equivalent to [q,p]* N = ih. With this 
quantization we must compute exp^ Ar (—*iLt//i) = exp * N (—iu>taa/h). The d- 
Schrodinger equation © actually is a simple first order partial differential equa¬ 
tion which shows, in particular, that the solution is a function of ad only. Writ¬ 
ing ad = s and denoting the solution by F(s,t ), the d-Schrodinger equation 
becomes 

ihd t F(s , t) = lusF(s , t ) + ujhd s F(s , t). 

The required solution, which must have the value 1 at t = 0, is F"(s,t) = 
e - s / h exp(e~ lUJt s/h), so we have 

exp * N (—iLutaa/h) = e~ aa ^ h exp{e~ lut aa/H), 


where on the right one has the ordinary exponential. Expanding the expression 
on the right and writing aa/h = H/Hco = p, the coefficient of e~ lnwt becomes 
e~ IJ, p n /n\. Comparing with © we see that the allowable values for the energy 
are E n = nhw and the corresponding = e~^p n /n\ (where N indicates that 
the normal form is used). We will write the (and generally those which 
arise with any quantization) as functions of H with the latter viewed as the 
energy function on phase space. For any value of H the are just the terms 
in the classical Poisson distribution with mean p = H/hui. This probability 
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distribution is sometimes called the ‘law of rare events’ 2 : if /t is the mean 
number of events seen in unit time (or space) then the probability that in a 
given unit of time (or space) one will see exactly n events is 7ri N ^ = e _#i /x”/n!. 
For example, if misprints are rare and the average number on a page is /i then 
the probability that a page will contain exactly n misprints is e _M /it”/n! (but 
/i may vary with the author). Unlike a Gaussian distribution, which depends 
on two parameters, its mean and standard deviation (‘uncertainty’ in physical 
terms, square root of its variance), the Poisson has but one, its mean. The 
variance is identical with its mean and the standard deviation is the square root 
of its mean. For large values of the mean, the distribution resembles a Gaussian 
with mean fj, and standard deviation ^//7. 

The question is how to interpret the appearance of the Poisson distribution 
here, bearing in mind that we have somehow the ‘wrong’ quantization (or at 
least not that in textbooks, since the lowest allowable energy is precisely zero, 
not hu/2). We will see that this quantization essentially presumes that we know 
the mean energy of the oscillator precisely, something which is not physically 
possible. In that impossible case it seems to say that if the mean energy is 
H then the probability that the oscillator will be observed in a state with the 
quantum number n is 7= e _A1 /x”/n!. But note that even though we have 
the ‘wrong’ minimal energy, the spectrum here is simply shifted by hio/2 from 
the textbook case, so the differences between allowable energy levels, which 
determine the spectrum, coincide with the usual. Nevertheless, the normal 
form of quantization is excluded by our selection principle. For with it all true 
probability distributions of energy are observable including a delta function, 
which has zero uncertainty, while the quantum uncertainty is always positive. 
(In fact, we will see that when the normal form is viewed as a limit the initial 
energy distribution is a delta function.) 

It is easy to verify that the 7r ™ are orthogonal idempotents in the *n 
multiplication, summing to 1 and having a fixed common integral over phase 
space: Direct computation shows that the integral of each is 2nh independent 
of n. That = 1 follows simply from setting t = 0. Finally, one 

way to see that the 7 r „ ^ are orthogonal idempotents is to observe that al¬ 
though *jv is a non-commutative multiplication, then^^ and exp* w (— icutaa/h) 
are all functions only of the single element ad and hence all commute. Com¬ 
paring the expansions of the two sides of the equation exp„ w (— iutad/h) *n 
ex P * N (~i <jJ taa/h) = eyLp* N (—2iutad/h) will show that 7 is idempotent. De¬ 
noting it for the moment by e, we clearly have for any idempotent e that 
e(l — e) =0, and that 1 — e is again idempotent. Proceeding by induction 
will show that the 7T« are all mutually orthogonal idempotents. 

While the normal form of deformation has been excluded by our selection 
principle, it already raises an interesting algebraic question equally meaningful 

2 Curiously, it arose neither from gambling nor physics but from Simeon-Denis Poisson’s 
study of the French judicial process in his “Recherches sur la probability des jugements en 
matiere criminelle et matiere civile”, 1837 



for all c-equivalent deformations. For the moment, let A denote the polynomial 
ring C [q,p] and let An denote the algebra to which we have deformed it. Note 
that as long as coefficients are restricted to polynomials in h it is meaningful 
to let h —> 0 in order to recover the original algebra; this gives the (only) 
correct statement of the correspondence principle. The resulting algebra is 
essentially the first Weyl algebra, a simple algebra (i.e., one without proper two- 
sided ideals) whose cohomology with coefficients in itself vanishes in all positive 
dimensions, in particular in dimension 2, and which is therefore rigid. With 
coefficients extended to C[[l//i]] [h] recovery of the original algebra by letting h —> 
0 is no longer possible. We have an algebra in which the identity has decomposed 
into a direct sum of infinitely many orthogonal idempotents e\, ... and whose 

precise structure we no longer know. The theorem that there is no degeneracy 
in one dimension suggests that each etAhe-i has dimension 1 and that the same 
is probably true for all eiAh^j- The simplest conjecture concerning structure 
would be that in each e t Ah.ej we can choose an element with eu = e ? ; and 
e-ijS-jk = 6ik, and that the algebra consists of linear combinations of these (but 
what beside the finite ones may be allowed is not clear). Again, the same 
question arises for all c-equivalent deformations of C [q,p\. 


3 Quantizations with A ^ 0 

Following the prescription in § 1, we now set T\ = e xhD D and define 

f *x 9 = T~\T x f * N T\ g) = f e *((i-A)££-A££> g . 

There are now several approaches to computing exp (—iHt/K). One can use 
the first equality above to get expressions for the quantities which we will now 
denote by , where 7r= 7 ; this will give them all in the form of power 
series. The second approach, which we adopt, is to use only the value obtained 
for 7 Tq A -* from the first method and then to adapt the procedure in 0 Appendix] 
to compute exp* (— iHt/K). The third approach is to note, as in the preceding 
section, that exp * (—iHt/h) will be a function only of aa; denoting this again 
by s and the result by F(s,t) one can solve the partial differential equation 
(now of second order in s) which F satisfies. However, we will see that one 
can also effectively solve the differential equation in closed form once closed 
expressions for the 7 have been obtained from the second method by using 
the generating function for the Laguerre polynomials. Since T\cia = aa + HA, 
with the first approach we have 

exp * x (-iHt/h) = 2)\" 1 exp* w {—iu>tT\{aa/h)) 

= e~ xhdadn exjp ifN {—iu>t{aa/h+ A)) 

= e- xiut / n e- xdada e- a ~ a/h e^ {e~ iuit aa/h) 

= e~ XiuJt ^ h e~ xdadii exp((e~ iut - 1 )aa/h). 


9 


A simple computation using the definition of the Laguerre polynomials then 
shows that 


exp , x (-iHt/h) = e~ Xiut J2(- X ) k ( e ~ iWt ~ l) k U(aa/\h) 

k—0 

where L*, is the fcth Laguerre polynomial. We may now replace aa/h on the 
right with H/tuo since the computations on the right now all take place in the 
undeformed algebra. Expanding the right side then gives 

OO OO / 7 \ 

exp ^(-iHt/h) = ^(-l)"^AM ? jL fc (L7/AMe“ i(n+A)a ' t - 

n—0 k =0 ^ 

The spectrum has thus been shifted, the allowable values of the energy are now 
E = (n + \)Tku, and we have 


= (-irjr\ n+k ( n t%n + k(H/\hv), 

k—0 x " ' 

which can be negative for some H. Using the generating function for the La¬ 
guerre polynomials, 


OO 

7u —1exp (rx _ ) = 5Z xA: ( -1 ) fcLfe ^)> 

i + x i + x f—; 

the special case n = 0 gives 


k—0 


^ = 


k—0 


( H ) 

1 ( 


VA Hu) 

1-A exp t 

(l — A )hujJ 


We can now adapt the procedure in [H] Appendix] to get from this the closed 
form for all the ir ^, namely 


7ri A ^ = 




Letting A —> 0 recovers original Poisson distribution. Comparing the two ex¬ 
pressions for gives the following identity involving Laguerre polynomials: 


(-i) n E A 


n-\-k 


k—0 


n + k 
k 

1 

~AVl- 


Ln+fc(^) 

-A \ n 




(1 - A). 


(9) 


(Setting A = 1/2 and z = 0 this asserts, for example, that J] fc (l/2) n+fc = 2 
independent of n, something easily verified directly.) Multiplying the left side 
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by e~ z / x lj m (z/\) = e -Z/,A L m ((l — A)z/A(l — A)) and integrating, the orthogo¬ 
nality relations of the Laguerre polynomials together with an obvious change of 
variables gives 


L m ((l - X)z)L n (z)e Z dz = 


if ™>n 
0 otherwise 


Replacing A by 1 — A and comparing the coefficients of the powers of A on the 
two sides gives the coefficients in the Fourier-Laguerre expansion of z k : 


z k L n (z)e 2 dz = 


= f (-!)”(*)*:! if k>n 
0 if k < n 


( 10 ) 


This fundamental result, which will show the existence of enough observable 
distributions, can also be derived in an elementary way, since it is just the 
formula for the change of basis from the Laguerre polynomials to the powers 
of x in the the vector space of polynomials in x: Let A be the infinite lower 
triangular matrix with rows and columns indexed by 0,1, 2,... and (i,j) entry 
equal to (—1 ) J ’Q), and D be the infinite diagonal matrix with diagonal entries 
1/n!, n = 0,1,2,... . Letting X be the infinite column vector (1, x, x 2 ,x 3 ,... )*, 
the nth entry in the vector L = ADX is just the Laguerre polynomial L n (x). 
To write X in terms of L it is sufficient therefore to invert AD , the only problem 
being the inversion of A. However, A is equal to its own inverse, for writing 
X r = (i - (i - x)Y = £(-!)* (iK 1 - %y gives 



1 if j = r 
0 otherwise 


Therefore X = D 1 AL, which is precisely what m asserts. 

While k is an integer in m~ one can deduce more generally that 


z p L n (z)e~ 


(-i r 


r(p + i ) 2 

nir(p — n + 1) 


With 01, the generating function for the Laguerre polynomials gives the follow¬ 
ing closed form for the *a exponential of — iHt/H , 


exp * (—iHt/h) = 


-iXujt 


■exp 


'2H e -iwt _ i 
1 - A + Xe~ iujt 




1 - A + Xe~ iujt 

at A = 1/2 one recovers the known formula for the Groenewold-Moyal case, 

1 


ex P*i /2 (~ iHt / h ) = 


cos(wf/2) 


exp (^L tan M/2))- 


The only values of A that need to be considered are 0 < A < 1/2. Setting t = 0 
shows that ^2 7T n'' > = 1 for any value of A; that the 7ri A ' 1 are mutually orthogonal 
idempotents follows exactly as in the case A = 0. We could, of course also deduce 
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this from the fact that the tt ^ are the transforms of the tt ^ by T \, but writing, 
as before, p, = H/Hoj we also have the duality relation 

TTi A) (/*) 7Tm _A) O) dn = 6 n ,ra- 
The Groenewold-Moyal case (A = 1/2) is self-dual. 



4 Negative probabilities and basic observable dis¬ 
tributions 


Returning to the question of negative probabilities, unlike the tt^\ we can 
not view the as giving an ordinary probability distribution over the energy 
values H since for all n > 0 and any positive A, will be negative for some 
positive value of H. (Moreover, we conjecture that for any fixed A / 0,1 and 
arbitrary H > 0 there must be an n > 0 such that < 0; this is easily seen to 
be true for sufficiently small H since ttq^ (0) = 1/(1 — A) > 0 and 7r "^ = 1 

for all H. In fact, while the foregoing sum is absolutely convergent for A < 1/2 
it is only conditionally convergent at the Groenewold-Moyal limit A = 1/2, and 
the convergence there is very slow.) 

Negative probabilities can not be dismissed as fiction. In the present case, 
accepting them at face value gives the correct expected value for the energy: 
differentiating the basic equation with respect to time and setting t = 0 
gives 

H = J2*eE, 

independent of the form of deformation. (Differentiating twice will give the 
quantum second moment rather than that of the original distribution, since the 
left side at t = 0 will be H * H .) In what remains we again write p for the 
dimensionless quantity H/hio when it is viewed as an ordinary scalar. 

As remarked at the beginning, one view of negative probabilities is that 
while they can not be observed directly we can observe positive distributions 
derived from them and thereby indirectly conclude their existence. Suppose 
that in our observation of the harmonic oscillator we have a true probability 
distribution p(p) for the energy, that is, one which is non-negative, defined for 
p > 0, and has / 0 °° p(fj,)dp = 1. With this, the probability of the oscillator 

being observed at the nth energy level becomes irl?\p,p) = / 0 °° n^\p)p(p)dp. 
For some probability distributions p these will all be non-negative; such p have 
been called observable distributions. The existence of sufficiently many is given 
immediately by us, for writing 


pi A V) 

one has from m that 
/*00 

/ n£\p)p[ X) (p)dp = 

Jo 


-J-) k e 
A fc! A ' 


Hl 

T 


(n) A"(l ~ A) fc ~ 
0 
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if fc > n 
0 if k < n 


The pare thus observable distributions. The coefficients on the right are 
the Fourier-Laguerre coefficients of the distribution (relative to the deformation 
with parameter A). We will call the p^ basic observable distributions since suit¬ 
able (in fact unique) linear combinations of them (necessarily involving negative 
coefficients) give all distributions with but a single non-zero Fourier-Laguerre 
coefficient (i.e., we can recover those given by the individual irh^). The basic 
observable distributions are the extreme elements of the convex cone of observ¬ 
able distributions and they span that cone. Note that the basic distributions 
have only a finite number of non-zero Fourier-Laguerre coefficients. 


5 The uncertainty inequality 

It is a classic computation that with the distribution p^ the mean or expected 
value of p, namely J pp^ (p) dp, is just A (k + 1). This is the mean value of the 
energy with the basic observable distribution p^ A . (Note that its minimum at 
k = 0 is, as expected, just A.) The second moment of the distribution is 

J p 2 p^\p)dp = X 2 (k + l)(k + 2 ). 

It follows that the variance is A 2 (fc+1), the square root of which is the standard 
deviation or ‘uncertainty’. The other evaluation of the expected energy in p^ 
necessarily gives the same result, for as observed earlier the energy calculation 
is independent of the quantization (and uses the negative probabilities). In fact, 
we have 

(n +\)ni l x \p)p[ x] {p)dp = J2( n +x )[ )A"(1- X) k ~ n = (k + l)X. 

Note that if we try to keep the mean of the distribution constant while letting A 
tend to zero, then the deformation becomes the normal one and the distribution 
becomes a Dirac delta supported at the mean energy. This is what was meant 
earlier by saying that the normal form of quantization assumes that the energy 
is precisely known. 

We can now apply our selection principle. With the basic observable distri¬ 
bution p^ one has 


OO r oo fc / 7 \ 

/ ( n + x ) 2 ^n ) (fj-)p i k\^)dp = ^2(n +\) 2 ( jA"(l — \) k ~ n 
n =0 n=0 A / 

= ( k 2 + k + 1)A 2 + kX. 

The variance now is fcA(l — A). By our selection principle, which asserts, in par¬ 
ticular, that quantization should not increase uncertainty, this quantum variance 
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must be smaller than the previous distribution variance. One therefore has the 
inequality 

fcA(l - A) < (k + 1)A 2 . 

We may not only assume that A is strictly positive, but as observed earlier, that 
0 < A < 1/2, so this implies that 


A > 


k 

2k+1' 


This being true for all non-negative k it follows that we must have A = 1/2, 
leaving the Groenewold-Moyal form as the only one consistent with the selection 
principle. Similar arguments may apply more generally to select one quantiza¬ 
tion from a family of cohomologically equivalent ones. 

Finally, fixing A at 1/2, note that the variance of the distribution p'^ 2 ' > is 
( k + l)/4 while its quantum variance is k/ 4, so the difference in variances is 
always just 1/4. The difference in uncertainties, however, is (Vk + 1 — Vk)/2 
which tends to 0 as the mean energy increases. This does not hold for any A 
strictly between 0 and 1/2, so for the simple harmonic oscillator it would be 
sufficient to require of a quantization procedure that as energy tends to infinitely 
(and quantization becomes unnoticeable) the difference in uncertainties tends 
to zero, but this might not be a sufficiently strong selection principle in general. 
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